Excitation of large-scale inertial waves in a rotating inhomogeneous turbulence 
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A mechanism of excitation of the large-scale inertial waves in a rotating inhomogeneous turbulence 
due to an excitation of a large-scale instability is found. This instability is caused by a combined 
effect of the inhomogeneity of the turbulence and the uniform mean rotation. The source of the large- 
scale instability is the energy of the small-scale turbulence. We determined the range of parameters 
at which the large-scale instability occurs, the growth rate of the instability and the frequency of 
the generated large-scale inertial waves. 
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I. INTRODUCTION 

The study of rotating flows is of interest for a wide 
range of problems, ranging from engineering (e.g., turbo- 
machinery), astrophysics (galactic and accretion discs) to 
geophysics (oceans, the atmosphere of the Earthy ra,seous 
planets) and weather predictions (see, e.g., 1, 2, 31). In- 
ertial waves arise in rotating flows and are observed in the 
atmosphere of the Earth and in laboratory rotating flows. 
In turbulent rotating flows inertial waves are damped due 
to a high turbulent viscosity. Thus, excitation of coher- 
ent and undamped inertial waves by turbulence seems 
not to be eflective. However, large-scale inertial waves 
are observed in turbulent rotating flows. A mechanism 
of excitation of the large-scale coherent inertial waves in 
turbulence is not well understood. 

Inertial waves are related with generation of large-scale 
vorticity. Generation of a large-scale vorticity in a heli- 
cal turbulence due to hydrodynamical a effect was sug- 
gested in This effect is associated with the aW 
term in the equation for the mean vorticity, where W are 
the perturbations of the mean vorticity and a is deter- 
mined by the hydrodynamical helicity of turbulent flow. 
A nonzero hydrodynamical helicity is caused, e.g., by a 
combined effect of an uniform rotation and inhomogene- 
ity of turbulence. 

Formation of large-scale vortices in a turbulent rotat- 
ing flows was studied experimentally and in numerical 
simulations (see, e.g., B IS E 111 IH El ) ■ Forma- 
tion of large-scale coherent structures (e.g., large-scale 
cyclonic and anticyclonic vortices) in a small-scale tur- 
bulence is one of the characteristic features of rotating 
turbulence (see, e.g., 0). A number of mechanisms have 
been proposed to describe generation of a mean flow by a 
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small-scale rotating turbulence, e.g., the effect of angular 
momentum mixing jl4j and vorticity expulsion 15]. The 
first experimental demonstration wherein it was shown 
that the divergence of the Reynolds stresses can gener- 
ate an organized mean circulation was described in jS] . 

There is a certain similarity between mean rotation 
and a mean velocity shear. Generation of a mean vortic- 
ity in a nonhelical homogeneous incompressible turbulent 
fiow with an imposed mean velocity shear due to an exci- 
tation of a large-scale instability was studied in [l&j . This 
instability is caused by a combined effect of the large- 
scale shear motions ("skew-induced" deflection of equi- 
librium mean vorticity) and "Reynolds stress- induced" 
generation of perturbations of the mean vorticity. This 
instability and the dynamics of the mean vorticity are 
associated with the Prandtl's turbulent secondary flows 
(see, e.g., US HI ES 111)- However, a turbulence with 
an imposed mean velocity shear and a uniformly rotating 
turbulence are different. In particular, the mean vorticity 
is generated by a homogeneous nonhelical sheared turbu- 
lence El- tlic other hand, the mean vorticity can- 
not be generated by a homogeneous uniformly rotating 
nonhelical turbulence (see below). The main difference 
between these two flows is that the mean velocity shear 
produces work in a turbulent flow, while a uniform rota- 
tion does not produce work in a homogeneous turbulent 
fiow. 

There are other interesting problems related with the 
inertial waves including, e.g., the effect of inertial waves 
on the onset of convection and on the turbulence dy- 
namics. In particular, the onset of convection in the 
form of inertial waves in a rotating fluid sphere were 
studied in 21]. On the other hand, the modiflcation of 
turbulence dynamics by rotation is due to the presence 
of small-scale inertial waves in rotating flows (see, e.g.. 

The main goal of this paper is to study large-scale 
structures formed in a rotating inhomogeneous turbu- 
lence. In particular, we investigate the excitation of 
large-scale inertial waves. These structures are associ- 
ated with a generation of a large-scale vorticity due to the 
excitation of the large-scale instability in an uniformly ro- 
tating inhomogeneous turbulence. The excitation of the 



2 



mean vorticity in this system requires an inhomogeneity 
of turbulence. 

This paper is organized as follows. In Section II we for- 
mulated the governing equations, the assumptions and 
the procedure of the derivation. In Section III the ef- 
fective force was determined, which allowed us to derive 
the mean-filed equations and to study the excitation of 
large-scale inertial waves in Section IV. The large-scale 
instability was investigated in Section IV analytically for 
a weakly inhomogeneous turbulence and numerically for 
an arbitrary inhomogeneous turbulence. Conclusions and 
applications of the obtained results are discussed in Sec- 
tion V. In Appendixes A, B and C the detailed derivation 
of the effective force is performed. 

II. THE GOVERNING EQUATIONS 



f jr VV ~ ^ ~ 

— = Vx(2Uxfi + :r(u)-t/VxW), (7) 

where ^(U) — —'^jifij — ffj^) is the effective force, 
fij — (uiUj) and f^j"^ is the second moment of the ve- 
locity field in a background turbulence (with a zero gra- 
dient of the mean velocity). Thus, the mean fields U 
and W represent deviations from the equilibrium solu- 
tion U*^"' = and W*^"' = 0. This equilibrium solution 
is a steady state solution of Eqs. Q and l@J. Note that 
the characteristic times and spatial scales of small-scale 
fluctuations of velocity and vorticity u and w are much 
smaller than that of the mean fields U and W. 

In order to obtain a closed system of equations in the 
next Section we derived an equation for the effective force 



The system of equations for the evolution of the veloc- 
ity V and vorticity W = Vxv reads: 



d_ 

dt 



V • V 



dW 



dt 



VP 

v = + 2 vx ft + lyAv + Fst , (1) 

P 

= Vx (vxW + 2vxJ7-i/VxW) , (2) 



where v is the fluid velocity with V • v = 0, P is 
the fluid pressure, Fgt is an external stirring force with 
a zero mean value, i7 is a constant angular velocity 
and ly is the kinematic viscosity. Equation (jSJ follows 
from the Navier-Stokes equation We use a mean 

field approach whereby the velocity, pressure and vortic- 
ity are separated into the mean and fiuctuating parts: 
v = U-|-u, P = P-|-p and W = W -I- w, the fluctuating 
parts have zero mean values, U = (v), P = (P) and 
W = (W). Averaging Eqs. (Q) and Q over an ensem- 
ble of fluctuations we obtain the equations for the mean 
velocity U and mean vorticity W : 



VP 

U = + 2Uxrj + :F + j/AU , (3) 

p 



dW 



Vx (UxVV^ + 2Uxrj+ (uxw) - z/VxW) , (4) 



where J-i 



Note that the effect of turbulence 



on the mean vorticity is determined by the Reynolds 

stresses (u,;u,) because 



(uxw)j = -Vj{uiUj) + 2^i(" 



(5) 



Consider a steady state solution of Eqs. © and I^J in 
the form: tj^**) = and W^'') — 0. In order to study a 
stability of this equilibrium we consider perturbations of 
the mean velocity, U, and the mean vorticity, W. The 
linearized equations for the small perturbations of the 
mean velocity and the mean vorticity are given by 



'dt 



VP 



(6) 



III. THE EFFECTIVE FORCE 

In this section we derive an equation for the ef- 
fective force The mean velocity gradients can af- 
fect turbulence. The reason is that additional essen- 
tially non-isotropic velocity fluctuations can be gener- 
ated by tangling of the mean-velocity gradients with 
the Kolmogorov-type turbulence. The source of energy 
of this "tangling turbulence" is the energy of the Kol- 
mogorov turbulence. The tangling turbulence was intro- 
duced by Wheelon 26] and Batchelor et al. 27] for a 
passive scalar and by Golitsyn (2^ and Moffatt 29] for a 
passive vector (magnetic field). Anisotropic fluctuations 
of a passive scalar (e.g., the number density of particles 
or temperature) are generated by tangling of gradients 
of the mean passive scalar fleld with a random velocity 
field. Similarly, anisotropic magnetic fiuctuations are ex- 
cited by tangling of the mean magnetic field with the 
velocity fluctuations. The Reynolds stresses in a turbu- 
lent flow with mean velocity gradients is another exam- 
ple of a tangling turbulence. Indeed, they are strongly 
anisotropic in the presence of mean velocity gradients and 
have a steeper spectrum (oc k~'^^^ ) th an a Kolmogorov 
turbulence (see, e.g., |3^|^|3^ 113,13]). The anisotropic 
velocity fluctuations of tangling turbulence were studied 
first by Lumley [s^] . 

To derive an equation for the effective force ^ we use 
equation for fluctuations u(f, r) which is obtained by sub- 
tracting Eq. for the mean field from Eq. Q for the 
total field: 



-2uxrj 



u 



N 



where 



U'^ = ((u • V)u) - (u • V)u -I- vAu 



(8) 



(9) 



We consider a turbulent flow with large Reynolds num- 
bers (Re = IqUq/v ^ 1), where uq is the characteris- 
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tic velocity in the maximum scale Iq of turbulent mo- 
tions. We assume that there is a separation of scales, 
i.e., the maximum scale of turbulent motions Iq is much 
smaller then the characteristic scale of inhomogeneities 
of the mean fields. Using Eq. (jS]) we derived equa- 
tion for the second moment of turbulent velocity field 
/y(k,R) = /(ui(k + K/2)uj(-k + K/2))exp(iK-R) dK: 



9/»j-(k,R) 
dt 



fijinn J mn 



fmn ^ij /) 



(N) 



(10) 



(see Appendix A), where Gijmn = ^ijmn(U) + Nijmn{^), the third moments / 



Equation H14() for the deviations of the second moments 
fij — ff^'' in k-space contains the deviations of the third 

moments fl^^ — fi^'^^ and a problem of closing the equa- 
tions for the higher moments arises. Various approximate 
methods have been proposed for the solution of problems 
of this type (see, e.g., [35,, 36, 37]). The simplest proce- 
dure is the T approximation which was widely used for 
study of different problems of turbulent transport (see, 
e.g., 35, 38, 39, 40, 41J). One of the simplest procedures, 
that allows us to express the deviations of the terms with 



(AT) 



fij^'^"^ in k-space in terms of 



9k- A A 



^im^jq^np 



^pU, , (11) 



N,,. 



2 ^qkpq ( 



'imp ^nj 



-jmp 



(12) 



and R and K correspond to the large scales, and r 
and k to the small scales (see Appendix A), Sij is the 
Kronecker tensor, kij — kikj/k'^, and V = d/dH, 
F,,(k,R) = (i?',(k,R)w,(-k,R)) + (M,(k,R)F,(-k,R)), 
F(k,R,t) = -kx(kxF,t(k,R))/fc2 and flp(k,R) are 
the terms which are related with the third moments ap- 
pearing due to the nonlinear terms. The third moments 

terms f-^^ are defined as 

)(k,R) = (F,„(ki);7„^(ki)«,(k2)) 

+ (,i,(ki)P,„(k2)J7„^(k2)) , 

where (k) is the Fourier transform of determined 
by Eq. ®, ki = k + K/2, kg = -k + K/2 and Py (k) ^ 

^ij kij . 

Equation H10|) is written in a frame moving with a lo- 
cal velocity U of the mean flow. In Eq. H10|l for the sec- 
ond moments of the turbulent velocity field we neglected 
small terms ~ O(V^) , where the terms with V ^ 0{L~^) 
contain the large-scale spatial derivatives. These terms 
are of the order of {Iq/LY, where the maximum scale of 
turbulent motions Iq is much smaller than the vertical 
size of the turbulent region L. 

Equation l|10|) for the background turbulence (with a 
zero gradient of the mean fluid velocity ViC/j = and 
r2 = 0) reads 



5/g\k,R) 
dt 



F^, + f, 



(N,Q) 



(13) 



where the superscript (0) corresponds to the background 
turbulence, and we assumed that the tensor i^,j (k, R), 
which is determined by a stirring force, is independent 
of the mean velocity gradients and of a constant mean 
angular velocity. Equation for the deviations fij — /^^ 
from the background turbulence is given by 



(0) 



dt 



ijmn T Jij J ij 



(14) 



that for the second moments fi^ 



reads 



[N] 



[N,0) 



fij ft 



(0) 



(15) 



where T{k) is the scale dependent correlation time of the 
turbulent velocity field. Here we assumed that the time 
T(fc) is independent of the mean velocity gradients (for a 
weak mean velocity shear). We considered also the case 
of slow rotation rate. In this case a modification of the 
correlation time of fully developed turbulence by slow 
rotation is small. This allows us to suggest that Eq. H15() 
is valid for a slow rotation rate. 

The T-approximation is different from Eddy Damped 
Quasi Normal Markowian (EDQNM) approximation. A 
principle difference between these two approaches is as 
follows (see ©El)- The EDQNM closures do not re- 
lax to the equilibrium, and this procedure does not de- 
scribe properly the motions in the equilibrium state. In 
EDQNM approximation, there is no dynamically deter- 
mined relaxation time, and no slightly perturbed steady 
state can be approached |35|. In the r-approximation, 
the relaxation time for small departures from equilib- 
rium is determined by the random motions in the equi- 
librium state, but not by the departure from equilib- 
rium ^3^. Analysis performed in ,35] showed that the 
T-approximation describes the relaxation to the equilib- 
rium state (the background turbulence) more accurately 
than the EDQNM approach. 

Note that we applied the r-approximation l|15|) only 
to study the deviations from the background turbulence 
which are caused by the spatial derivatives of the mean 
velocity and a uniform rotation. The background turbu- 
lence is assumed to be known. Here we use the following 
model for the background isotropic and weakly inhomo- 
geneous turbulence: 



/f (k,R) 



Snk^ 



I 

2k 



2 (fc^V, 



-fcjVi) 



,2\(0) 



(16) 



(see e.g., |43|), where Py(k) — Sij — hj, Sij is the 
Kronecker tensor and kij = kikj/k'^, T{k) — 2Tof(fc), 
£{k) = -df{k)/dk, f{k) = (fc//co)^~«, 1 < g < 3 is the 
exponent of the kinetic energy spectrum (e.g., q — 5/3 
for Kolmogorov spectrum), fcg = l/lo- 



4 



The mean velocity gradient ViU causes generation of 
anisotropic velocity fluctuations (tangling turbulence). 
Equations (|14|I -H16 |I allow to determine the second mo- 
ment fij (R) = / fij (k, R) dk: 



(see Appendix B), where 



(0) 



6 

To(u2)(o)/6 and the 



(17) 



tensors My and 5^ are determined in Appendix B. 
The definition of the function yields (u2)(")(R) = 
Gvj^^IVj/To. Since we assumed that tq is independent of 
R, the spatial profile of the function 1^^°^ [e.g., given by 
Eq. H2()|l in Sect. IV-B] determines the spatial profile of 
(u^)'^'^^. Equation 1171) allows to determine the effective 
force .F. = -V,[/„ (R)-/{;'^(R)] : 

(18) 

where A = (V/2)//2 ^ (Vi/(,"))/z.(,o). 

Note that when A = const, the effective force does 
not have the term cx aW, where a describes the hy- 
drodynamic a-effect. The hydrodynamic a-effect was 
introduced in the equation for the mean vorticity (see, 
e-g-, ^, IE 01) J similarly to the a-effect in the equa- 
tion for the evolution of the mean magnetic field (see, 
e.g., 01) • The reason for the absence of the aW term 
in is as follows. Let us suggest the opposite, i.e., 
that ^ oc aW = aVxU. Since the effective force 
^^ = -V,[/i,(R) - /f°'(R)], we obtain 



-aeijkUk 



(19) 



Here we used the identity Wi — Sijk^ jUk and we took 
into account that when A — const, the hydrodynamic a 
is constant. Note also that in our paper we considered 
incompressible velocity field. The condition H19|) is in 
contradiction with the Galilean invariance, because the 
Reynolds stresses in the considered case may depend on 
the gradient of the mean velocity field rather than on 
the mean velocity itself. When A is not constant, the 
effective force J- can have the term c>c aW. However, 
this effect is not in the scope of our paper (e.g., this case 
cannot be described in the framework of the gradient 
approximation) . 



IV. THE LARGE-SCALE INSTABILITY IN AN 
INHOMOGENEOUS TURBULENCE 

For simplicity we consider the case when the tur- 
bulence is inhomogeneous along the rotation axis, i.e., 
A = A(z)e^, Q, = VlBz- After calculating [Vx(VxU)]2 
from Eq. © and Wz from Eq. Q we arrive at the fol- 
lowing equations written in nondimensional form 




0.1 0.2 0.3 0.4 0.5 Q.x„ 



FIG. 1: The rotation rate dependence of the functions I3{Q.to) 
(solid) and /3^^(f2ro) (dashed). 
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FIG. 2: The rotation rate dependence of the functions 
v^iyiTo) (solid), u^I^tq) (dashed-dotted), ri^(Q.T{)) (dashed) 
and ri^(Q.Ta) (dotted). 



dt 



-f i/„ AAi V, -f A^V^];/, , 
= (G-/3AA_l)?7^ + KV^ + i^^AV, 



+ 77„,AJVF, , 
where A = A^ + V^, = d/dz, 

G = 2a,V, +/3[AV2-f AV^] 



(20) 



(21) 



AAV, 



I3{n,z 
v^{Vl, z 

v^j(yL,z 

r]^{n,z 



(z)(../8)[i?5(c^)/2-7^6(u;)] , 
u!^Hz){iu/8)D,{u;) , 
u!^Hz)[D,{u;)/2 + D2{u;)] , 
r^,{n,z) + 2v'^^\z)D2{Lo) , 
v^^\z)D^{uj) , 
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is given by 



Here uj = StqO, a* = fiL^/z^* > 1, and we used Eq. (fTH|l . 
The functions Dk{uj) are determined in Appendix B. 

In Eqs. (|20l) and 1)21(1 we use the following dimension- 
less variables: length is measured in units of L, time in 
units of L^/i/* , the function A(z) is measured in the units 
of i^^, the function v^\z) is measured in the units of 
^'*, the perturbations of velocity Uz and vorticity Wz 
are measured in units of J7* and U<,/ L, respectively. The 
functions I3{VLtq), I3^{^Itq), v^{VLtq), i^^{nTo), ?7^^.(f^To) 
and r]^{^lTo) are shown in FIGS. 1-2. All these func- 
tions shown in FIGS. 1-2 are normalized by i>'^\z), e.g, 

^^{flTo) = ^)/'^!^^H^)j similarly for other func- 

tions. 



IPinro) A\K- K(r!To) + V^inro)] 



(25) 



The maximum growth rate of perturbations, 7,„ — 
[P{flTo) A]'^ /&[i^^{nTo) + r]„{nTo)], is attained at if = 

= IPinro) A\/2[i^^{nTo) + Vu^i^T-o)]- This case cor- 
responds to a* ^ P"^ {Qto) / Av^ {Qtq) . The large-scale 

hydrodynamic helicity of the flow is U • (VxU) = 

4ifj.|U.|2sgn(/3) ^0. 



B. Numerical results 



A. The weakly inhomogeneous turbulence 

Assume that functions and A{z) vary slowly 

with z in comparison with the variations of the mean 
velocity Uz{z) and mean vorticity Wz{z). Let us seek 
for the solution of Eqs. H2U|) and (|21|l in the form oc 
exp(7i — iK • R). Let us first consider perturbations with 
the wave numbers <^ K^. Since V • U = 0, the ve- 
locity components Uz <C |U^|. Thus the growth rate of 
the inertial waves with the frequency 



uJw = -sgn(/3 A) ■ 



2 a» AT, 



K 



is given by 



In 



\f^{nTo)AKz\-iy^{nTo)K^ 



(22) 



(23) 



where J — Jw + i^un the wave number K is mea- 
sured in the units of L^^ and 7 is measured in the 
units of v*/L^. The maximum growth rate of the in- 
ertial waves, 7m = [/3(l^ro) A]^/4i^y (firo), is attained at 
K = Km = |/3(Oto) A|/2z^y (r^To). For a very small rota- 
tion rate, i.e., for to = S^Itq <^ 1, the turbulent viscosity 
ly^inro) ~ (g+3)/5 and Pinro) « (32/15)17ro, where the 
parameter q is the exponent of the kinetic energy spec- 
trum of the background isotropic and weakly inhomo- 
geneous turbulence (e.g., g = 5/3 for Kolmogorov spec- 
trum), and this parameter varies in the range 1 < q < 3. 
Note that the inertial waves are helical, i.e., the large- 
scale hydrodynamic helicity of the motions in the inertial 
waves is tj • (VxU) = 2\lJ^\'^Kz ^ 0. This instability is 
caused by a combined effect of the inhomogeneity of the 
turbulence and the uniform mean rotation [see the first 
term in Eq. 

Now we consider the opposite case, i.e., the per- 
turbations with the wave numbers ^ K^. Since 
V • U = 0, the velocity components Uz ^ |Uj_|. When 
K'^ |/3(ri) A|/4a* < < K±, the growth rate of per- 
turbations with the frequency 



u-u, = sgn(/3 A) 



2 a* Kz 
K 



In this Section we take into account the inhomogeneity 
of the functions v':^\z) and A(z). We introduce a new 
variable V = ISUz and consider an eigenvalue problem for 
a system of Eqs. (|20f) and 1)21(1 . We seek for a solution of 
Eqs. (|2DJ) and (EB in the form oc '^{z) exp{'jt)jQ{Kj^ r), 
where Jo{x) is the Besscl function of the first kind. After 
the substitution of this solution into Eqs. I(20|l and l(21|l we 
obtain the system of the ordinary differential equations 
which is solved numerically. 

We used the cylindrical geometry (z,r, 0) with z-axis 
along the rotation axis and consider the axisymmetric 
solution (i.e., there are no derivatives with respect to 
the polar angle (f). The turbulence is inhomogeneous 
along the rotation axis. We use the periodic bound- 
ary conditions in z direction for Eqs. ((2()(l and H21|l. i.e., 
C7,(z = 0,r) = f7,(z = L,r), C/^(z = 0, r) = f7^(z = L, r), 
il'^iz = 0,r) = U'^iz = L,r), U'^'iz = 0, r) = C/f (z = 
L,r), Wz{z = 0,r) = Wz{z = L,r) and W^{z = 0,r) = 
W!,{z = L,r), where U' — ^ zU ■ We also use the con- 
dition Ur{z,r = 0) = Ur{z,r = B) = 0, where R is the 
radius of the turbulent region. 

We have chosen the vertical profile of the function 
v^^z) in the following form 



i/(°)(z) = l-C<il-exp 



C 



-2A 



z 1 
L " 2 



(26) 



1 - 



1 - exp(-A2/2) ' 



(24) 



with two values of the parameter v'^'' = 0.1 and 0.3; and 
two values of the parameter Aq = 12 and 30. The ver- 
tical profile of the turbulent viscosity i^^*''(z) is shown 
in FIG. 3. The maximum of turbulence intensity is lo- 
cated at z = L/2. The form of the chosen spatial profile 
of the function ly^^ (z) is simple enough and universal. 
It allows us to vary the size of the region occupied by 
turbulence (by changing the parameter Aq) and the dif- 
ference in the level of the turbulence between the center 
and boundary of the region (by changing the parameter 
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z/L 



0.2 



0.4 



0.6 
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1 



FIG. 3: The vertical profile of the turbulent viscosity u!^'' (z) 
for Ao = 30 and ^''^ = 0.1 (solid); Ao = 30 and 4'' = 0-3 
(dashed); Ao = 12 and i/^*"' =0.1 (dashed-dotted). 
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1^^^), i.e., it allows us to change the inhomogeneity of 
the turbulence. The numerical solution of Eqs. H20() 
and l|21(l was performed also for other spatial profiles of 
the fimction i'!:^'>{z). However, the final results do not 
depend strongly on the details in the spatial profile of 
the function i^^*'-'(z). Note also that the chosen spatial 
profile of the function i''^^ (z) can mimic the distribution 
of turbulence in galactic and accretion discs (see, e.g.. 

The sufficient condition for the excitation of the insta- 
bility is 7„ > 0. The range of parameters {L/Iq, Qtq) for 
which the large-scale instability occurs is shown in FIG. 4 
for different values of the parameters /i, iy^^{z) and Aq. 
Here /i — L/Lr, L is the vertical size of the whole re- 
gion, Lr is a radius from the center of the structure at 
which the energy of the radial velocity perturbations 
is maximum. Note that the maximum radial (horizon- 
tal) size R of the whole region is of the order of ^ ALr- 
The decrease of the parameter ^ causes increase of the 
range of the large-scale instability. On the other hand, 
the increase of the size of the highly intense turbulent 
region (i.e., decrease of the parameter Aq) results in the 
increase of the range of the instability. 

The rotation rate dependencies of the growth rate 
"twTj^ of the large-scale instability and the frequency uJwTo 
of the generated waves due to the large-scale instability 
are shown in FIGS. 5-6, where r^^ = L'^/v*. There is a 
threshold in the rotation rate for the large-scale instabil- 
ity: r2*ro « 0.025, and when 17 > f^.^, the instability is 
excited. The instability threshold in the parameter L is 
L > IOIq. 

Note that the characteristic time (~ 2Tr/jw) of the 
growth of perturbations of the mean fields U and W is 
by 5 orders of magnitudes larger than the turbulent cor- 
relation time tq. The period of oscillations T = 2n/uJio 
of inertial waves is at least 10 times larger than the tur- 



FIG. 4: The range of parameters {L/Io,^Itq) for which the 
large-scale instability occurs (7^, > 0) for (a). Ao = 12, 
4"' = 0.1; (b). Ao = 30, = 0.1; (c). Aq = 30, 4''> = 
0.3; and different values of the parameter fi: fi — 0.1 (dotted), 
fi — 0.5 (solid), /J. — 1 (dashed), /j. — 2 (dashed-dotted). 



bulcnt correlation time tq. The minimum value of the 
period of rotation = 27r/f2 is at least 20 times larger 
than the turbulent correlation time tq. All spatial scales: 
the vertical size of the turbulent region, L, and the ver- 
tical size of the highly intense turbulence (~ L/Aq) are 
much larger than the maximum scale of turbulent mo- 
tions Iq (e.g., L/Aq is at least 10 times larger than the 
maximum scale of turbulent motions lo)- This implies 
that there is indeed a separation of scales as we assumed 
in the derivations. Note also that the range of validity of 
the obtained results is tq ^ T for a statistically station- 
ary background turbulence. In particular, we assumed 
that the characteristic time of evolution of the second 
moments is much smaller than the turbulent correlation 
time. The asymptotic formulas (|22|I - H25() are in agree- 
ment with the obtained numerical results. 

It must be noted that turbulent flow with an im- 
posed mean linear velocity shear and uniformly rotat- 
ing flows are essentially different. In particular, in a 
turbulent flow with an imposed mean linear velocity 
shear there are no waves similar to the inertial waves 
which exist in a uniformly rotating flows. The reason 
is that any shear motions have a nonzero symmetric 
part {dU)ij of the gradient of the mean velocity, where 
{dU)ij — (yiUj + VjC7i)/2. In addition, the difference 
between these two flows is that the mean velocity shear 
produces work in a turbulent flow, while a uniform ro- 
tation does not produce work in homogeneous turbulent 
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FIG. 5: The rotation rate (Qto) dependencies of (a) the 
growth rate ^wT,-, of the large-scale instability and (b) the 
frequency cjwTq of the generated waves due to the large-scale 
instability for Ao = 12, f^''' = 0.1, = 0.1 and different 
values of the parameter L/Iq: L/Iq = 50 (solid), L/lo = 
100 (dashed) and L/lo = 500 (dashed-dotted). Here Tj^ = 



FIG. 6: The rotation rate (Qto) dependencies of (a) the 
growth rate ^wT^ of the large-scale instability and (b) the 
frequency uiwTo of the generated waves due to the large-scale 
instability for Aq = 30, f!^^ = 0.1, ^ = 0.1 and different 
values of the parameter I///o: Z///o = 100 (solid), L/Iq =500 
(dashed) and L/lo = 1000 (dashed-dotted). 



flow. 



V. CONCLUSIONS 

We studied formation of large-scale structures in a ro- 
tating inhomogeneous nonhelical turbulence. We found 
a mechanism for the excitation of the large-scale iner- 
tial waves which is associated with a generation of a 
large-scale vorticity due to the excitation of the large- 
scale instability in a uniformly rotating inhomogeneous 
turbulence. It was shown that the mean vorticity can- 
not be generated by a homogeneous uniformly rotating 
nonhelical turbulence. The excitation of the mean vor- 
ticity in this flow requires also an inhomogeneity of tur- 
bulence. Therefore, the large-scale instability is caused 
by a combined effect of the inhomogeneity of the turbu- 
lence and the uniform mean rotation. The source of the 
large-scale instability is the energy of the small-scale tur- 
bulence. The rotation and inhomogeneity of turbulence 
provide a mechanism for transport of energy from turbu- 
lence to large-scale motions. We determined the range of 
parameters at which the large-scale instability occurs. 

Some of the results obtained in this study, e.g., the 
expression for the effective force in a homogeneous tur- 



bulence, are in compliance with the previous studies of 
rotating turbulence 47] [see Appendix B, Eq. ljB13p ]. It 
is plausible to suggest that the results of recent experi- 
ments |48| can be explained by the large-scale instability 
discussed in this paper. 
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APPENDIX A: DERIVATION OF EQ. STHi 

In order to derive Eq. H10(l we use a mean field ap- 
proach, i.e., a correlation function is written as follows 

(ui(x)itj(y)) = y"(uj(ki)uj(k2)) exp[i(krx 
-fks-y)] dki dk2 
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= J Ji-j (k, R) exp(ik-r) dk , 

/„(k,R) = y"(u,(k + K/2)7i,(-k + K/2)) 

X exp(iK-R) dK 

(see, e.g., E^); where R and K correspond to the 
large scales, and r and k to the small scales, i.e., R — 
(x + y)/2, r = x-y, K = ki + ka, k = (ki - 
k2)/2. This implies that we assumed that there exists a 
separation of scales, i.e., the maximum scale of turbulent 
motions is much smaller then the characteristic scale 
of inhomogeneities of the mean fields. 
Now we calculate 

9/,j-(ki,k2) _ ,p .cKi(ki) 

= (nn(kl) Wj(k2)) 

+(^.(kOP,„(k2)^^^) , (Al) 



dt 



where we multiplied equation of motion ||SJ| rewritten in 
k-space by -Py (k) = 5ij — kij in order to exclude the 
pressure term from the equation of motion, Sij is the Kro- 
necker tensor and fcy = kikj/k'^. This yields the equation 
for /ij (k, R) [see Eq. ((TUI)]. For the derivation of Eq. ifTUI) 
we used the following equation 

iki j fij{k - iq, K - Q)C/p(Q) exp(iK-R) dKdQ 



« f'dU 



4, ,|)(V.V.C/,). 



dk. 



(A2) 



To derive Eq. (|A2|) we multiply the equation 
V • u = 0; written in k-space for Mi(ki — Q), by 
Mj(k2)f/p(Q) exp(iK-R), and integrate over K and Q, 
and average over the ensemble of velocity fluctuations. 
Here ki = k + K/2 and k2 = -k + K/2. This yields 



1 



1, 



i [k, + -K, - Q,j {u,{k + 2^ - Q)wj 
X UpiQ) exp (iK-R) dK dQ = . 



1 



K)> 



(A3) 



Next, we introduce new variables: ki = k + K/2 — Q, 
k2 = -k + K/2 and k = (ki - k2)/2 = k- Q/2, K = 
ki + k2 = K — Q. This allows us to rewrite Eq. (jA3|) in 
the form 

i (^h + \k, - Uj (k - ^Q, K - Q)C/p(Q) 
X exp (iK-R) dK dQ = . (A4) 
Since |Q| <C |k| we can use the Taylor series expansion 
/,,(k - Q/2, K - Q) ~ /,,(k, K - Q) 



i a/,,-(k,K-Q) 

2 dk. 



(A5) 



We also use the following identities: 



/,,(k,R)C/p(R)]K = / /.j(k,K-Q);7p(Q)dQ, 



Vp[/y (k,R)C/p(R)] = j zifp[/y(k,R)[/p(R)]K 
X exp (iK-R) dK , (A6) 

where [/^ (k, R)L/p(R)]K denotes a Fourier transforma- 
tion. Therefore, Eqs. HA4p - HA6p yield Eq. IIA2p . 

APPENDIX B: THE REYNOLDS STRESSES 

In this Appendix we derive equation for the Reynolds 
stresses using Eq. (I14II . We assume that the character- 
istic time of variation of the second moment (k, R) is 
substantially larger than the correlation time T(fc) for all 
turbulence scales. Thus in a steady-state Eq. (|14|l reads 

[L{n) - T{k) /(u)] (/ - /(")) = T{k) [N{n) 

+/(u)] /(o) , (Bl) 

where we used Eq. (|15|l . Hereafter we use the fol- 
lowing notations: / = /y (k,R), /(^) = 4'^^(k,R), 
/(O) = /^^(k,R), /(U)/ = /,,„„(U)/„„(k,R) and 
N{n)f = 7V,,>„„(0)/,„„(k,R), and L(f2) = i,y™„(n) = 
5im53n - T{k) N^jrnni^)- Multiplying Eq. ||BT1) by the 
inverse operator L^^ijl) yields 

[E - r(fc)L-i(n)/(U)](/-/(o)) = -[^-L-i(f2) 
-r(fc)L-i(0)/(U)]/(°) , (B2) 

where E = SimSjn and we used an identity 

E - L-Hn) = -T(fc) L-\n) N{n) . 

The latter identity follows from the definition: 
L^^{fl)L{fl) = E. The inverse operator L~^{fl) is 
given by 



" — 1 — 1 

(r2) = ijj,,„„(r2) = -[Bl SimSjn + B2 kijmn 



^"-^5 ^ipni^jqnkpq ~\~ Bq {Sipmkjpji -\- Sjpnkipni^] 7 (-^^) 

where Bi = 1 + x(2V'), B2 = Bi + 2 - 4x(V'), 
B3 = 2V'x(2V'), Bi = 2xW - Bl, B5 ^ 2- Bi, 
Be = 2^[x(V') - X(2V')], x{x) = 1/(1 + and ^ = 

2T(fc)(k • n)/k. 

Multiplying Eq. ljB2p by the operator 
E + T{k) L~^{i},) I (IJ) yields the second moment 
/^/,,(k,R): 



/ « [L-'{n) + T{k)L-\n)i{ij)L-\n)]f^'\ 



(B4) 
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where we neglected terms which are of the order of 
0(|VU|2). Since Lr.^JO) F„„(k) = /',,(k), Eq. m 
reads 



/ « /(")+r(fc)L-i(f2)/(U)/ 



(0) 



(B5) 



The first term in Eq. (jBSp describes the background tur- 
bulence. The second term in Eq. (|B5|I determines effects 
of both, rotation and mean gradients of the velocity per- 
turbations on the turbulence. The integration in k-space 
yields the second moment fij (R) — J fij (k, R) dk which 
is determined by Eq. (|17|l . where we used the notation 

/ = fij — fij ~ fij 

given by 



and the tensors Mij and Sij are 



+DiiLo){uj-V)iu,-i])u;,j , (B6) 
S,, = D^{lo) + Dq{uj) i?y , (B7) 



where 



Qij = {bJiVj + u)jVi){G} ■ U) 

+ {Cj-V){u:,U,+ufii) , (B8) 
T^j ^ (tZ;xV).,(a>xU)j + ((2)xV)j(i2>xU), , (B9) 



(BIO) 



Ri 



[w,,(tDx V) 



,(i2)xV)j]((2) • U) 



+ (w • V)p,((2)xU)j +c2;j(a;xU),] , (Bll) 

{dU)ij ~ (S/iUj + \7 jUi)/2, ujij — iliiLOj, u)i — ^li/VL, and 

D,{u) = {2[A«H-4^)(0) + (q + 2)cf^(0) 
+cWH]+A«H}/4, 



^ [2dl\u)~A^^>{u)]/S, 



(2), 



uj = Sto^I. The functions Am\w) and Cm\uj) are de- 
termined in Appendix C. Equation (|18|l for the effective 



(2)/ 



force ^ = — Vj/ij(R) can be rewritten in the form 

= T^"^ + A/,, + 5,, ) A, , (B12) 

where = to(u2)(")/6 = ^oVStq, A = (V;§)//2 = 

(Vi/(°))/i^(,"), and the effective force in a homoge- 

neous turbulence reads 



{H) 



AO) 



+ {D2 + D^)[{CJ-VfU, 
1 



l2)^ A(t2> • U)] 



De (oD- V)[(a)xV),(cD-U) 
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-(a>-W)w, + (tD- V)(a>xU)i] 



(B13) 



where we used the identity 



^ijk^lmn — ^il^jjn^kn "t" ^in^jl^km ^ini^jn^kl 

— SinSjmSkl — SilSjnSkTn — SimSjlSkn ■ 

Equation (|B13|) for the effective force J-", in a homoge- 
neous turbulence coincides in the form with that obtained 
in 47] using symmetry arguments. However, the symme- 
try arguments cannot allow to determine the coefficients 
in Eq. (|BT3|) . 



APPENDIX C: THE IDENTITIES USED FOR 
THE INTEGRATION IN k-SPACE 

To integrate over the angles in k-space we used the 
following identities: 



Jy(a) = 



kij sin 9 
1 + a cos^ 9 



d9d(p = AiSij + A2UJij , (CI) 



J 



J ijrnn (^) 



ijmn 



(a) 



^ijrnn sill 

1 + a cos^ 9 



ijmn 



(1 



= J, 



ijnin 



(a) 



da 



dO dip = 
(a) 



^ijmn sin ^ 

6 + a cos^ 6 



dd dip 



b=l 



(C2) 



(C3) 



where = (2)^(2;^, ujijmn 
^2 = C2 -f 7C3, and 



I ^ ^ 27r 
^i(a) — 
a 



(a+l) 



i^ij^mn, Ai — 5Ci -|- C3, 

arctan(Y^) ^ 



A2ia) 



2tt 
a 



arctan(V^) ' 
(a -I- 3j -= 6 
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,,arctan(->/a) 5a 



(Sa^ + 30a + 35) 



arctan(y^) 



(53(a) 



55a 
'~3~ 

TT 

2^ 



35 



(a^ + 6a + 5) 



arctan(y^) 13a 



y/^ 3 
In the case of a ^ 1 these functions are given by 

Ai{a) ~ (47r/3)[l - (l/5)a] , ^(a) - -(87r/15)a , 
Ci{a) ~ (47r/15)[l- (l/7)a] , C'2(a) ~ (327r/315)a2, 
(73(a) (87r/105)a. 

In the case of a » 1 these functions are given by 

Ai{a) ~ TrVVa, ^2(0) Tr^Va , 

Ci{a) ^ TT /4^/a - 47r/3a , C2 (a) ~ 37rV4\/a , 
(73(a) n'^/Ay/a + Sir/Sa . 

Now we calculate the following functions 



,,P+i 



(a;) = (6/7ra;^'+^) / y^Ck{y')dy 



The integration yields: 

A^P{uj) = 12 



arctan(u;) ( ^ 1 



+J_[l_ln(l+a;2)] 



A«(a;) = -12 



arctan(w) A 3 



:[3-21n(l + w2-)] 



arctan(a;) 



3-4.-i 



Af)(a,) = 
Af{u:) = 
Cf\uj) = (3/2) 



arctan(u;) 



1 \ 3 2 ; 
1 + — - — + — 



arctanfw) 1 \ 7 6 ^, , 
— 1 + — - — + — 



arctan(u;) ( ^ 1 



13 
3^ 



A^i\Lo)-7A^l\u)+i^&^\Lo) , 
A^^\^)-5C^f\u.), 



where S{u)) = Jg [arctan(y)/2/] dy. In the case of a; <C 1 
these functions are given by 



AW(a;) 


~ 4(1- 


la;^),AW(a;)~ 


4 2 
-5" ' 






A^2),Af(c.)~ 


16 2 

a; , 

25 ' 


C«(-) 


4 

- 5^^- 


lc.2),c(^\c.)^ 


0(u;4) , 




4 

<~ a; 

35 


^Cf)(a;)~A(l 


35 ^ 




~ 0(a;*) 


,Cf(c.)~-J|a;^ 
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